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Within the framework of liquid crystal flows, the Qian & Sheng (QS) model [22] for Q —tensor 
dynamics is compared to the Volovik & Kats (VK) theory [26] of biaxial nematics by using 
Hamilton’s variational principle. Under the assumption of rotational dynamics for the Q— tensor, 
the variational principles underling the two theories are equivalent and the conservative VK 
theory emerges as a specialization of the QS model. Also, after presenting a micropolar variant 
of the VK model, Rayleigh dissipation is included in the treatment. Finally, the treatment is 
extended to account for nontrivial eigenvalue dynamics in the VK model and this is done by 
considering the effect of scaling factors in the evolution of the Q— tensor. 


1 Introduction 

Following the celebrated Ericksen-Leslie (EL) model 
[7, 21] of uniaxial liquid crystals, several dynamical 
theories have been formulated over the decades. These 
theories often differ by the choice of order parameter 
(depending on the phase under consideration) and even 
different director formulations are available in the sim¬ 
plest case of uniaxial molecules. For example, while 
the EL theory defines the director as an unsigned 
(position dependent) unit vector in physical space, the 
Harvard theory [12] defines the same order parameter as 
a differential one-form or vector field, as elucidated in 
Volovik’s work [27]. This difference leads to intrinsically 
distinct transformation properties of the director field, 
which are reflected in substantial differences between 
the EL and the Harvard theories. 

In the description of other nematic phases, such as 
the biaxial phase, de Gennes [3] showed the convenience 
of identifying the order parameter with a traceless 
symmetric tensor, the Q— tensor. This has the advan¬ 
tage of incorporating features of different liquid crystal 
phases (e.g. uniaxial and biaxial) by simply writing 
different expressions of the same tensor order parameter 
in terms of the director fields. Again, different theories 
of Q —tensor dynamics are available. Other than by the 
expression of the free energy or dissipation function, 
some of these theories again differ by intrinsic trans¬ 
formation properties of the Q— tensor. For example, 
in [21 ] the alignment tensor is a (degenerate) covariant 
tensor field, while in Qian & Sheng [22] the same tensor 
is simply a (position dependent) 3x3 matrix. As 
before, this leads to substantial differences between the 
theories. In other situations, alignment tensor theories 
differ in whether or not inertial effects are retained in the 
rotational dynamics. When these effects are neglected, 
dissipation dominates to drive relaxation dynamics. For 
example, inertial rotation terms are neglected in the 


Q —tensor dynamics of Volovik & Kats [26], while the 
theory of Beris & Edwards [1] neglects inertial effects 
to formulate Q —tensor dynamics as a gradient flow on 
the space of symmetric covariant tensor fields. More 
generally, dissipation can be introduced by either the 
use of symmetric brackets [1] (which are added to the 
Poisson bracket governing inertial effects) or the use of 
Rayleigh dissipation in Hamilton’s variational principle 
[24, 6], 

In recent years, different approaches to director 
dynamics were shown to be equivalent when the director 
is defined as a unit vector. Indeed, while it is well known 
how the EL theory can also be formulated in terms 
of the molecular angular velocity (see Volovik’s works 
[5, 26]), new developments [15, 13] have shown how 
several liquid crystal theories are actually special cases 
of Eringen’s theory of micropolar fluids [8, 9], Based on 
the variational approach, variants of these theories are 
constructed depending on how one expresses the free 
energy and the dissipation function. As a result, these 
studies allowed the identification of the hydrodynamic 
helicity invariant for the conservative limit of several 
liquid crystal dynamical models [16]. 

In this paper, we compare the Q —tensor dynamical 
theories of Qian & Sheng (QS) [22] and Volovik & Kats 
(VK) [26]. More particularly, upon restoring inertial 
effects in the VK theory (see also [17]), we show that 
its dissipationless limit arises as a specialization of the 
QS model, under the common assumption of rotational 
dynamics for biaxial liquid crystal phases, that is [18] 

Q = RQ q R t (1) 

(where R(x,t) is a rotation matrix). This evolution 
naturally arises from de Gennes’definition (see e.g. [ 3]) 

Q = (u u - f(u)d 2 u, ( 2 ) 

where / is the statistical distribution of the molecular 
unit vector u (pointing along the long molecular axis) 
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on the unit sphere S 2 . We shall show that using the 
above evolution in the variational principle for the QS 
model returns VK theory augmented to allow for inertial 
effects. More particularly, upon following Edwards’ 
approach in [6], most of this paper focuses on the 
use of Hamilton’s principle to formulate conservative 
dynamics, and this picture is eventually extended by the 
insertion of Rayleigh dissipation [24] . After discussing 
the variational principle for the conservative limit of the 
QS model for liquid crystal textures, Section 2 makes 
use of the evolution (1) to recover the VK theory in the 
absence of fluid flow. In addition, Eringen’s wryness 
tensor is introduced as an auxiliary variable to obtain 
a micropolar variant of the VK model. Then, Section 3 
extends the treatment to flowing liquid crystals. Section 
4 deals with Rayleigh dissipation upon making use of 
the evolution (1) eliminating the need for the constraint 
TrQ = 0 appearing in QS theory. Finally, Section 5 
extends the treatment to consider nontrivial eigenvalue 
dynamics for the Q —tensor and this is done upon 
invoking the existence of a scaling factor, by mimicking 
the properties of microstretch fluids [10, 11, 14]. 

2 Conservative texture dynam¬ 
ics: QS and VK theories 

In this section, we show that the QS model specializes 
to the VK theory, under the assumption of rotational 
evolution (1). Upon considering the special case of a 
flowless nematic texture in the absence of dissipation, 
we shall derive the QS and VK theories from Hamilton’s 
variational principle, along the lines of Edwards [6] . In 
particular, we shall make use of the rotational symme¬ 
try by applying standard Euler-Poincare techniques in 
geometric mechanics, following the works of Holm [19] 
and Gay-Balmaz & Ratiu [14]. 

Upon denoting by J the moment of inertia density, 
the conservative QS model for Q —tensor dynamics in a 
nematic texture reads 

Jd 2 t Q-h= A 0 1 (3) 

where h is the molecular field 


field such that h = h T . The QS equation (3) is evidently 
an Euler-Lagrange equation on the space of symmetric 
matrices, arising from the constrained variational prin¬ 
ciple 

sj y (j IIQII 2 - HQ, VQ) + Ao Tr q) d 3 x dt = 0 , 

(4) 

where the last term is the constraint enforcing the 
traceless condition. Recently, a similar approach was 
followed for biaxial phases in [18], where constrained 
variational principles were used to prescribe the values 
of the first three rotational invariants in the sequence 
{TrQ™}. On the other hand, in the same work, it is 
recognized how biaxial molecules require the Q— tensor 
to undergo purely rotational evolution of the type (1). 
In again the same paper, the authors emphasize how 
the change of position of a molecule in space can be 
described by the angular velocity matrix, which in our 
notation reads 

v = [d t R)R~ x 

and whose corresponding angular velocity vector is Vi = 
— eikjVjk/2- It is the purpose of this paper to introduce 
this variable explicitly in the QS model and show how 
the resulting equations coincide with the VK model in 
[26], when the latter is augmented to consider inertial 
effects arising from J ^ 0. An immediate verification of 
this statement is obtained upon noticing that (1) implies 
dtQ = [y, Q], which is then replaced in (3): then, taking 
the commutator of the latter equation with Q yields the 
conservative VK equation [26] Jdt[\y,Q],Q] = [/i,Q]. 
The remainder of this paper shows that the reason 
for the appearance of these equations lies in the fact 
that the QS and VK models share the same Hamilton’s 
variational principle, under the ansatz (1). Therefore, 
their underlying minimization problems are equivalent. 

In order to introduce the variable v in the variational 
framework, we apply Euler-Poincare theory [20, 19, 14] 
by replacing the evolution relation (1) in the action 
principle (4). This operation yields the Euler-Poincare 
variational principle 

sj:j (j II [y, Q}\\ 2 d 3 x - HQ , VQ)) d 3 * dt = 0, 


h = 


dT_ 

dQ 


+ div 



J- is the Landau-de Gennes free energy 


where the constraint has now been dropped appro¬ 
priately and the rotational invariants Tr Q n no longer 
contribute to the Landau-de-Gennes free energy. Then, 
upon computing the Euler-Poincare variations 


•F — 2 i^Qij “I” Qij t k + RlQij,jQik,k) 

47t _ 

V Ll^ijkQil^jQkl [dQijQjkQki T "f{Qij) 

with a, L\ 1 L/ 2 , /3,7, P phenomenological constants, and 
Ao is a Lagrange multiplier enforcing Tr Q = 0, so that 
taking the trace of (3) yields 

Ao = ^Tr(Jd 2 Q-h) =~\^h. 

In the original paper [22] , an extra Lagrange multiplier 
enforces Q = Q T , although here we drop the corre¬ 
sponding term by suitably symmetrising the molecular 


Su = d t rj + [rj, v ], SQ = [rj, Q\ , (5) 

with rj = 5RR- 1 , one is led to the Euler-Poincare 
equations 

J^- t [[y, Q],Q] = [h, Q], ^ = [y, Q\ ■ ( 6 ) 

which coincide with the conservative limit of the VK 
model in [26], upon setting J = 0 to neglect inertial 
effects (and in the absence of fluid flow). These effects 
are also recovered in the VK model by retaining the 
corresponding terms in the Poisson bracket structures 
appearing in [26]. Notice that the traceless condition 
has never been imposed in this treatment because 
it naturally arises from a convenient choice of initial 
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condition Tr Q o = 0. This Euler-Poincare variational 
approach was recently followed in [2] to study textures 
of biaxial nematic phases in an external field. 

Perhaps not surprisingly, the above equations com¬ 
bine into 

\J'Q ~h,Q] =0, 

which is solved, for example, by 

Jd 2 t Q -h = \ n Q n , 


VK model of biaxial liquid crystals, we shall proceed to 
consider the more general case of flowing liquid crystals. 
In this case the rotational degrees of freedom are cou¬ 
pled to the relabelling properties that characterize fluid 
flows. 

3 QS and VK theories for liquid 
crystal flows 


for arbitrary time dependent functions \ n . Alterna¬ 
tively, if Q had not been considered a priori as obeying 
(1), then X n would correspond to a sequence of Lagrange 
multipliers prescribing the values of TrQ n+1 . As an 
example, in [18], the first three A„’s were considered, 
while the QS theory uses Ao only. Consequently, the 
VK model emerges as a more reliable alternative to the 
approach followed in [18], since the rotational dynamics 
is intrinsically encoded in the model. Indeed, the Euler- 
Poincare formulation of the VK model was recently 
adopted in [2], 

In conclusion, we have proved that the conservative 
VK theory emerges from the QS model, under the 
assumption (1) of purely rotational dynamics of biaxial 
molecules. This relation between the two theories arises 
from the fact that both are derived from exactly the 
same action principle (4). 

A micropolar theory of the above equations can also 
be obtained by noticing that [15] 

d z Q = [diRR-\Q\ + RfyQoR- 1 
= [d l RR- 1 ,Q\ + R[Q 0 ,i 0i ]R- 1 
= [Q, -diRR- 1 + RjoiR- 1 } 

=■ [Q,%\ 

where we have invoked the existence of an initial wry- 
ness tensor such that diQo = [Qo,7oi]- In Eringen’s 
micropolar theory [9], the time-dependent wryness ten¬ 
sor 

7i ; = ^ d l RR ~ 1 + 

identifies the amount of rotation under an infinitesimal 
displacement dx and thus determines the spatial rota¬ 
tional strain [20]. Upon computing 

$li = [V, 7i] - dv 


As in the previous section, we start by focusing on the 
following QS equations for liquid crystal dynamics: 

P D t u = -Vp - dt (^V,,) , (9) 

pJD\Q -h = A 0 1 (10) 

D t p + pdivu = 0 (11) 

where D t = dt + u ■ V is the ordinary convective 
derivative. Notice that the above equations are a slight 
generalization of the original QS equations in [22], as 
they allow for a compressible fluid flow and the pressure 
is expressed in terms of an internal energy function 
U ( p ) as p = p 2 U' — T. It is straightforward to observe 
that equations (9)-(10) arise from the following Eulerian 
action principle: 

s 1 y(>i 2 +>^ii 2 -^) 

— T(Q, VQ) + Ao TrQ^j d 3 xdt = 0 , (12) 

with variations 

5u = dfW + (u ■ V)w — w ■ Vtt 
Sp = — di v(pw) 

5Q = —(w ■ V)Q + 0 
6(D t Q) = -(w ■ V)D t Q + D t O 

for arbitrary w and 0 vanishing at the endpoints. 

Although the form of the variations above may look 
somewhat mysterious, they find a natural justification in 
terms of the relabelling properties of the action in (12). 
Here we shall follow an analogous treatment to that in 
[13]. In order to show how this works, we shall introduce 
the Lagrangian fluid path X(xo, t) and its corresponding 
Lagrange-to-Euler map 


and by writing F(Q, [Q,7]) = ^(Q, 5 /), one obtains the 
equations of motion 


d_ 

dt 

dQ 

dt 


= [u Q ], 





n 

. dQ 

+ 



d 

(7) 


dji 

dt 


= [£,7i] - diV. 


( 8 ) 


These equations represent the micropolar version of 
the VK model. In the gauge theory of defects, an 
inhomogeneous initial condition on the wryness tensor 
can be associated with the presence of disclinations in 
the texture [4, 20, 13]. 

Now that we have characterized how the QS and VK 
theories are related in the case of conservative texture 
dynamics, and have provided a micropolar variant of the 


p(x,f) = J p 0 (x 0 )<5(x-X(x 0 ,t))d 3 a;o =: X*p 0 , 

whose fundamental role is expressing the Eulerian den¬ 
sity p(x, t ) in terms of its (fixed) Lagrangian correspon¬ 
dent po( x o)- Then following Euler-Poincare theory we 
write the Eulerian velocity as u(x.,t) = X(X' 1 (x, t), t) 
and observe that the relabelling property of the action 
in (12) takes it into the form: 

jf/Qpo||X|| 2 + ^ollQf - Po^(X* Po o X) 

- P °, v - R(Q, Vg) + (X.Ao) Tr q) d 3 xdt (13) 

(X*po)°X J 

where we have introduced Q = Q o X and o denotes 
standard composition of functions. In turn, the above 
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action identifies the expression of a Lagrangian L po of 
the type 


similar fashion to that of the previous section, the 
Q —equations combine into 


L P 0 (X,X,Q,Q)=e(p,u,Q,D t Q), (14) 


\pJD 2 Q-h,Q} = 0, (18) 


where £(p,u,Q, D t Q) is the Lagrangian appearing in 
the Eulerian action principle (12), which has now been 
written in (13) in terms of purely Lagrangian variables. 
The relation (14) means that the Lagrangian description 
(arising from L po ) is equivalent to the Eulerian descrip¬ 
tion (arising from £) - in fluid mechanics, this is typically 
known as ‘relabelling symmetry’. 

At this point, without entering into the difficult 
question of computing the Euler-Lagrange equations 
for X and Q, we observe that the previous Eulerian 
variations arise naturally from the definitions above. 
Thus, one has 

Su = S(X o X -1 ), SQ = S(Q o X -1 ), 5p = 6{X,p 0 ) 

along with the relations w = (<5X) o X~ 3 , Q = (SQ) o 
X” 1 and (SfS)oX " 1 = D t Q. 

Now that we have unfolded the relabelling features 
of the QS model, it is easy to proceed analogously to the 
previous section by making the evolution ansatz (1) for 
biaxial liquid crystal flows. Indeed, with the definitions 
above, we write 

Q = TlQQRr 1 (15) 

(where lZ(x.o,t) is a rotation matrix) and introduce the 
angular frequency matrix u> = (dt'R-)'R.~ 1 . Then, the 
variations SQ and 8u> are expressed as in (5) (upon 
replacing 9 by uj and Q by Q). However, in order to 
obtain the corresponding Eulerian description, we define 


which is the restriction of equation (10) under the 
assumption (1) of rotational evolution. 

In conclusion, we have proved that the conservative 
VK theory of biaxial liquid crystals is a specialization 
of the QS model under the assumption (1) of purely 
rotational dynamics, and that this is due to their 
common action principle (12). 

A micropolar fluid version of the above VK equations 
can be obtained by following the same steps as in 
the previous section, upon defining Q := —d{\ZTZr 1 + 
Hit'll- 1 and 7 such that 


PI ■= / Po(x 0 )f(xo)<5(> 


X(x 0 , t)) d 3 a; 0 . 


The equations of motion can be written upon computing 

7 i + (u- V)7i + (diu k ) 7^ = [9,7i] - dp) 

5 % + (w ■ V)7i + ( diw k )% = [77,7i] - dpp 


and by writing J-(Q, [<3,7]) = ^(<9,7). Then, upon 
denoting p = p 2 U’ — T, the action principle (16) yields 


pD t Ui 


-d iV - ftTr 




D t Q = [9, Q\ 

(19) 


pJD t [[u,Q],Q] = 


D t % + ( diu k )% = 


' -TV 



Q) dQ 

+ 

-^,1 i 
_0^fi 


[9,7*] - dpv 



( 20 ) 


( 21 ) 


Q = {TZQ 0 n 1 )oX 1 = QoX 1 These equations represent the micropolar version of the 

9 = ( d t un _1 ) o X" 1 =uio X~ 3 . VK model for flowing liquid crystals. 


so that we find 

D t Q = [9, Q] 

SQ = [f?, Q\ - (w ■ V)Q 

8v = d t fj — (w ■ V)9 + (u ■ V)rj + [rf, 9] 

with rj = (STZTZ~ 1 )oX~ 1 arbitrary and vanishing at the 
endpoints. Then, Hamilton’s principle associated to the 
action (13) is transformed into its Eulerian formulation 
as 

6 1 7" + Q] II 2 

-pU{p) - D(Q,X7Q)J d 3 xdt = 0 (16) 

and by using the variational relations above, one is led 
to the Euler-Poincare equations (9), (11) and 


4 Rayleigh dissipation and rota¬ 
tional dynamics 

This section introduces dissipation in the previous liquid 
crystal models. Although dissipation can be introduced 
in the Hamiltonian framework by the use of symmetric 
brackets [ ], the present framework makes the use of the 
Rayleigh dissipation function £%(Q,Q), in which case 
the variational principle associated to the QS equation 
(3) for texture dynamics becomes 

sj t 2 J (j WQW 2 - HQ , VQ) + Ao Tr q) d 3 * dt 

-jCMi 40 )** (22) 

so that (3) changes to 


pJD t [[9, Q],Q] = [h, Q] , D t Q = [9, Q] (17) 

Again, these equations coincide with the conservative 
limit of the VK model in [26], upon setting J = 0 
in order to neglect inertial effects, and they do not 
involve Lagrange multipliers since the invariants Tr Q n 
are naturally preserved by the evolution of Q. In a 


Jd 2 Q — h = AqI — 


8&_ 
SQ ' 


Here we have used the notation for the functional 
derivative, whose definition reads 


< 5 F(w) 


-— • (5wd 3 a; 
Sw 
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On the other hand, under the rotational evolution 
(1) for biaxial nematics, one has &(Q, Q ) = r(y, Q ) 
with the variations (5), so that the relation S£% = Sr 
yields 

Sr St% 

~SV = [SQ ,Q . ' 

Eventually, as noticed already in [13], the Euler- 
Poincare variational principle with dissipation 


which then replaces the right hand side of (12) 

so that 

the QS equations (9)-(ll) become 



( dJF \ 

Sr 


pD,u- Vy> a( a% v<3«) 

“ Su’ 

(27) 

pJD t Q h-Xol S ( DtQ y 


(28) 

D t p + pdivu = 0 . 


(29) 


1 2 


^||[9,Q]|| 2 -.F(Q,VQ))d 3 ;rdt 


42 r / x r 

Tr ( — rf) d 3 z dt (23) 


yields the equations of motion 


r*| 

dt 1 


these combine into 


Sr 

SV 


dQ 

dt 



Sr 

.Q], 

SV 


Then. 

conser- 


Again, 



S2% ' 

JQ — h H-— Q 

SQ 


= 0, 


Alternatively, one may proceed by assuming the evolu¬ 
tion (15) in &p 0 to write 

/^ 0 (X,X, Q. Q) d 3 x 0 = J r(u,p,V, Q) d 3 i, 

with (25) and 

^(xo,t), 2(x 0 ,t) S(x - X(x 0 ,t))d 3 x 0 ■ 
SQ 


5@, 


P 0 


SX 


SX + Tr 




Sr 


= /( — •in + Trl—^ 77 ) d 3 a; 


Su 


Sr 


SV 


where the dot notation stands for partial time deriva¬ 
tive. The same approach applies to the micropolar 
theory (7)-(8), whose dissipative version is obtained by 
inserting the term —SrjSV into the right hand side of 

(7). 

In the case of flowing liquid crystals, dissipation can 
be introduced by mimicking the same steps as above. 
For a consistent theory, it is essential to start with the 
Lagrangian L po (X,X,2, 2) defined in (14) and write 
its dissipative action principle 


L Po (X, X, Q,Q)dt 


=0(lF' ^x+ K1H) Ad, 


where 2# Po = ^ Po (X,X, Q, Q). Then, upon assuming 
that depends on X only through p = X*po 

and using the same notation and definitions as in the 
previous section, we can write 

^ P0 (X,X,Q, Q) =r(u,p,Q,D t Q ) 


with 


Sr f S& 

-f—=J—^r-{xo,t)S(x-X(x 0 ,t))d 3 x 0 , (25) 

r 

S(m = J^t M S { x-X(xo,t) ) d 3x o. (26) 


and the dissipative version of the VK equations is 
obtained by inserting the term —Sr/SV into the right 
hand side of the first of (17), which is then accompa¬ 
nied by (27) and (29). After verifying that Sr/SV = 
[5r/S(D t Q),Q\, these equations combine into 


pD t u = -S/p - di 


dT 

dQij j 

Sr 


pJD/Q h+ 5{DtQ) ,Q 
D t p + pdivu = 0 . 


S/Qij 
= 0, 


Sr 
Su ’ 


At this point, the micropolar formulation of VK 
dynamics can also be presented in dissipative form 
by following precisely the same steps as above. This 
consists of adding the term —Sr/Su into the right hand 
side of the first in (19) and the term —Sr/SV into the 
right hand side of (20). 

It is important to emphasize that different Rayleigh 
functions produce different relaxation dynamics and 
therefore different physical results. Here, we shall not 
dwell upon the difficult question of what the correct 
expression of the Rayleigh function should be; see [24], 
Still, the results presented here remain valid indepen¬ 
dently of the particular choice of Rayleigh functional. 


5 Scaling and eigenvalue dynam¬ 
ics 


Eventually, 


8& P 0 

SX 


■ SX -|- Tr 


' SM pt 

V <52 


■<52 


Sr 


= I [ — w- 1- Tr 
Su 


d 3 x 

Sr 


S(D t Q) 


SQ d 3 * 


The evolution (1) is isospectral - that is, the eigenvalues 
of Q are preserved. However, more general cases are 
also allowed in liquid crystal dynamics. Then, the 
purely rotational motion can be extended in order to 
incorporate nontrivial eigenvalue dynamics. In the 
present geometric context this is achieved by allowing 
the evolution law to admit scaling. In light of this, 
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one choice is to let Q evolve under dilations as well as 
rotations, that is 

Q = X 2 RQ q R t , (30) 

where A(x, t) is a time-dependent positive function. 
For example, in the case of an initial uniaxial phase 
Qo = cr(nono — 1/3) (for some constant a), the above 
relation yields 


where a\ 2 may now acquire the meaning of a scalar 
order parameter. A more general situation can be 
studied by introducing two conformation tensors A and 
B such that Q = X 2 RAR T + ( 2 RBR T . However, this 
possibility is left for future work. 

The evolution law (30) comes out in a natural 
fashion by considering the congruence transformation 
Q = gQog T for some invertible matrix g(x,t). Indeed, 
as shown in Appendix A, a natural choice is a matrix 
g of the form g = A R. As an additional remark, we 
recall that such conformal rotations play a crucial role 
in Eringen’s description of microstretch fluids [10, 11], 
as enlightened in [14], This motivates us to investigate 
the interplay between scaling and rotational degrees of 
freedom. To this purpose, we simply extend the Euler- 
Poincare theory from previous sections to accommodate 
dilations. 

Upon using the evolution law (30), one follows the 
procedure outlined in Section 3 by defining 

Q = Q o X- 1 = (X 2 TZQ a 'Rr 1 ) o X“\ 
v = (d t nn- 1 ) ox- 1 , 

A=(A' 1 9 ( A)oX- 1 , 

so that the variational principle (16) is modified to 

Q p\M 2 + ^p\\W, Q]|| 2 + 2pJA 2 ||Q|| 2 

-pU(p) - R{Q, VQ)j d 3 xdt = 0, (31) 

with the variations 


Q = aX 2 


SQ = [% Q ] + 2SQ - (w ■ V)Q (32) 

Sv = d t rj — ( w ■ V)9 + ( u ■ X7)rj + [77, 9], (33) 

SA = d t E - (w • V)S + (u ■ V)A, (34) 

where 77, E and w are arbitrary and vanish at the end¬ 
points. Then, the resulting Euler-Poincare equations 
extend the VK model to incorporate scaling: 


pD t u = -Vp - di 


dR 


dQ 


ij ,1 


VQj 


D t p = —pdivit, 
(35) 


P JD t [[9, Q],Q] = [h, Q ], D t Q = [9, Q\ + 2AQ, 

(36) 

2pJ||Q|| 2 (AA + 2A 2 ) = pJ||[9, Q ]|| 2 + TV (hQ ), (37) 


where we observe that the first three equations remain 
unchanged, while the modified Q —tensor dynamics pro¬ 
duces an additional equation for the scaling quantity 


A (the latter quantifies the time rate of Tr Q n far all 
integers n). Also, notice that the dynamics above can 
be rewritten upon replacing (36) and (37) by 

[pJD 2 Q ~h,Q} = 0 , Tr ({pJD 2 Q - h)Q) = 0 , 


where the latter is obtained by combining (37) with 
(36). Note that the above dynamical system is also 
a special case of the QS model, since the former is 
obtained by making use of the ansatz (30) in the 
variational principle (12) underlying the QS equations 
of motion (9)-(ll). 

At this stage, inserting dissipation is straightforward 
as it can be done by simply applying the method 
outlined in the previous section. Therefore, we shall 
just state the final results, which are as follows: 


pD t u = — Vp - di 
pJD 2 Q — h + 


dR 

&Qij ,1 

Sr Q 


VQi, - 


S(D t Q) 


JpTr(QD 2 Q) = Tr( [h- 
D t p + pdivw = 0 . 


= 0 , 
Sr 


S{D t Q) 


Sr 

Su 



As before, the Rayleigh dissipation function is left arbi¬ 
trary to allow for different modeling options. Possible 
choices of Rayleigh function are presented in [21]. 


6 Conclusions and open ques¬ 
tions 

This paper has compared three apparently different 
approaches to Q —tensor dynamics in the theory of 
liquid crystals. Starting from simple texture dynamics, 
the VK theory for biaxial nematics was shown to emerge 
from the QS model under the preliminary assumption of 
conservative dynamics and by retaining inertial effects. 
In addition, a micropolar variant of the VK theory has 
been formulated by following Eringen’s definition of the 
wryness tensor [9]. Already at this stage the use of 
Euler-Poincare reduction theory is advantageous, and 
even more so in the case of flowing liquid crystals, which 
was the subject of Section 3. The relation between 
QS theory and VK dynamics was shown to hold in 
this more general case without essential modifications 
and the same holds for the corresponding micropolar 
variant. Finally, we applied Euler-Poincare reduction 
to the dissipative action principle, thereby showing how 
rotational dynamics ( 1 ) can also be taken into account 
for dissipative liquid crystal flows. Notice that all the 
results obtained in this paper also apply to the special 
case of incompressible flows. 

Another natural question was addressed regarding 
the possibility of Q —tensor flows that are more general 
than ( 1 ). Indeed, while an isospectral flow of Q is 
justified for biaxial molecules, other types of nematic 
molecules require more general evolutions which are 
then influenced by the terms Tr Q n in the Landau-de 
Gennes free energy. A possible strategy to tackle this 
question was pursued by making use of the conformal 
rotation group, so that (1) is replaced by (30). This 
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approach was followed by Eringen in his theory of 
microstretch fluids [10, 11], and its geometric features 
have been exploited in [14]. In the present paper, a 
full dynamical theory was developed based on the Qian- 
Sheng Lagrangian and dissipation was also included. 

The work in this paper follows in the same direc¬ 
tion as recent work [13], which showed how Eringen’s 
micropolar theory includes Ericksen-Leslie dynamics as 
a special case. It is hoped that the identification of 
relations between different liquid crystal theories will 
help provide a better understanding of the essential 
features captured in the various models. As mentioned 
in the introduction, an outstanding question concerns 
possible deep relations between the Harvard theory and 
EL dynamics. The answer to this question is among 
the most difficult, since the order parameter possesses 
essentially different transformation properties in the two 
models. 
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A Dilations, rotations and traces 

This appendix shows how conformal rotations may 
arise from the transformation properties of traceless 
symmetric matrices. In particular, we shall show the 
following 

Lemma A.l Let g G GL(n,R) be such that 

gQg T G sym 0 (3), VQ € sym 0 (3), 

where sym 0 (3) denotes the space of traceless symmetric 
matrices. Then, g is of the fom 

g = \R 

with A G R + and R G SO(3). 

Proof. Clearly the quantity gQg T , g G GL(3) is sym¬ 
metric. To preserve tracelessness one has ( g T g,Q ) = 0 
as an inner product on sym(3). Since this form is 
nondegenerate, g T g must be an element of sym(j-(3), 
the orthogonal complement to the traceless subspace of 
sym(3). Then noting the following: 

dim sym(3) = 6 
dim sym 0 (3) = 5 
dim sym^(3) = 1, 

inferred by counting, the only symmetric and tracefull 
1-dimensional space orthogonal to every element of 
sym 0 (3) is A 2 1 for A 2 G R, since g T g is positive definite. 
Therefore g = A R for R G SO(3). ■ 
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